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Abstract 

It is generally recognized that economical systems, and more in general complex 
systems, are characterized by power law distributions. Sometime, these distributions 
show a changing of the slope in the tail so that, more appropriately, they show a 
multi-power law behavior. We present a method to derive analytically a two-power 
law distribution starting from a single power law function recently obtained, in the 
frameworks of the generalized statistical mechanics based on the Sharma-Taneja- 
Mittal information measure. In order to test the method, we fit the cumulative 
distribution of personal income and gross domestic production of several countries, 
obtaining a good agreement for a wide range of data. 

Key words: Two-power law distribution, Sharma-Taneja-Mittal information 
measure, distribution of personal income and gross domestic production. 
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1 Introduction 



Free-scale behavior in the economical systems have been observed since 19th 
century, when Pareto noticed that the cumulative distribution of the personal 
income P(x) = f£°p{y) dy of several countries behaves like a power law func- 
tion. Afterwards, Gibrat clarified that such a power law behavior holds only 
for the high income region, whilst in the low-middle income region, which in- 
cludes almost the whole body of data, the curve is well fitted by a log-normal 
distribution. 

Actually, the problem concerning the real profile showed by the function P(x) 
in the whole range of the accessible data is still an open question. In particular, 
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it has been suggested pQ that deformed exponential functions derived recently 
in the field of the generalized statistical mechanics, can be fruitfully employed 
to modeling analytically the cumulative distribution P(x) for a wide range of 
the income values. 

Notwithstanding, the recent analysis based on a huge quantity of data nowa- 
days accessible, shown that sometime the crossover among the low-middle 
region (the log-normal region) and the high region in the upper tail of the 
distribution (the Pareto region, with a power law behavior P(x) ~ x~ s , where 
s is a positive constants quite generally 1 < s < 2), does not occur smoothly, 
giving origin to knee or ankle effects (see for instance [2]). Moreover, in some 
cases, it has been observed a deviation from the Pareto behavior in the high- 
est region, which can originate a new power law behavior P(x) ~ x~ s with a 
different slope s^s. 

The complicate profile in the shape of P{x) cannot be accounted for by a 
generalized exponential with a single power law behavior. This open the ques- 
tions: how can we describe the shape observed in P(x) with an analytically 
simple function? 

In the present contribution, we introduce a mechanism which permits to gen- 
erate multi-power law functions by employing deformed exponentials and log- 
arithms with a single power law asymptotic behavior. 

Notice that, two-power law behavior have been observed in various econom- 
ical systems like, for instance, in the cumulative distribution of the personal 
income [3], in the cumulative distribution of the land price [I] or in the returns 
of many market indexes [5]. 

On a general ground, two-power law behavior have been observed in different 
physics fields as well as in biological, geological and social sciences. Among 
the many, we quote the dielectric relaxation [6], the re-association in folder 
proteins [7], and others [8]. 

It is worthy to remark that there have been proposed different methods in 
literature [Tt^yTU] to produce generalized distributions with a double-power 
law behavior which differ from the one advanced in the following. 



2 Deformed logarithms and exponentials 

Generalized exponential functions S(x), interpolating between the standard 
exponential exp(x) for x < 1 and the power law x~ s for x ^> 1, arise naturally 
in the study of thermostatistic proprieties of complex systems which show free- 
scale feature. In [IT], it has been postulated a very general expression for the 
entropy of such a system 
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(in the unity of Boltzmann constant k-g = 1), where A(x) plays the role of 
a generalized logarithm, the inverse function of 8{x). By requiring that the 
distribution, derivable through a variational problem, assumes the form 



M 



v{x) = e\-Y,Pj 



\x 



(2.2) 



0=1 



which mimics the well-known Boltzmann-Gibbs distribution, the following 
functional equation has been obtained 



d 

dx 



xA(x)} = A A ( 
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Here, a and A are constants given by 
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The quantities (3j in Eq. ( 12. 2D play the role of Lagrange multipliers associated 
to the M constraints / |x| Mj 'p(x) dx = Oj which represent the //j-th momenta 
of x. Typically, the constants fij are integers (for instance, n\ = gives the 
normalization / p(x) dx = 0\, [12 = 1 is the mean value (x) = O2, and so on) 
but for sake of generality we assume fij e M. 

The most general solution of Eq. (12.31) . accounting for the boundary condi- 



tions A(l) = and (d/dx) A(x) 



x=l 



1, derived from certain physically and 



mathematically justified assumptions, is given by 



x 



X — X 

2k~ 



(2.5) 



which recover the standard logarithm in the (k, r) — > (0, 0) limit. 
By requiring that ln {t _ r} (x) is a continuous, monotonic, normalizable, con- 
cave and increasing function for x G (0, +00), we obtain the restrictions 
-|«| < r < \k\, if < \k\ < 1/2 and |k| - 1 < r < 1 - |k|, if 1/2 < \k\ < 1. 
Notwithstanding, for particular applications some of the above mathemati- 
cal requirements can be relaxed permitting less restrictive conditions for the 
deformation parameters. For instance, in certain practical situations one is 
welling with a normalization in a finite interval x G (0, x max ) [12] and we can 
discard the condition \k\ < 1. 

In the following, we require only that Eq. ( 12. 5ft be a monotonic function, so 
that its inverse function, the generalized exponential exp| K r \(x), certainly ex- 
ists. This is accomplished by requiring only that — \k\ < r < \k\. 
From Eq. (I2.5P we obtain that ln {Kr} (x) — > x r+ ' K '/|2fi;| for x — > +00 and 
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ln {K r y(x) — > — x r l K l/|2 k\ for x — > 0, whilst ln {K r} (x) — > (a; — 1) for | ii? — 1 1 <C 1. 
In the same way, we have that exp| K r \(x) — > |2 k a^VC^M) for a; — > ±oo, whilst 
exp| K r |(x) — > 1 + x for x — > 0. Thus, the deformed exponential expr K r y(x) 
interpolates with continuity between the standard exponential exp(x) ~ 1+x, 
for x — > 0, and the power law \x\~ s with slope s = — l/(r ± \k\), for x — > ±oo. 
Finally, accounting for the solution (12.51) . the entropy (12.11) assumes the form 

S .,r(P) = ~ P( X ) ln {K , r} (P( X )) dX , ( 2 - 6 ) 



which recovers, in the limit (k, r) — > (0, 0), the Shannon-Boltzmann-Gibbs 
entropy S = — J p(x) kip(x)dx. This entropic form, introduced previously 
in literature in [T3fl4|ll5] . is known as the Sharma-Taneja-Mittal information 
measure and has been applied recently in the formulation of a possible ther- 
mostatistics theory [TfiyTT] . 



3 Two-power law function 

Endowed with the deformed logarithm ln {K r} (x) and the deformed exponential 
exp { (x) we can construct the quantity 

ri CTl (x) = exp {Ki n} (ai ln {Ki ri} (x)) , (3.1) 

where o\ denotes the set of parameters <j\ = (ki, r%, a\), with a\ > 1. The 
function (13. ip is therefore employed in the following construction 

f(x) = U ai o exp ra} (x) = U ai (exp (-x)) . (3.2) 



We observe that, for a x = 1 expression (13.21) reduces to exp { ^ r (— x), for 
(«i, ri) — ► (0, 0) we obtain [exp { ^ r (— x)] ai , whilst for (/c 1; n) = (k 2 , r 2 ) we 
obtain exp r , (— Oi a?). 

Accounting for the asymptotic behavior of the deformed exponential and log- 
arithm we can distinguish three regions in the range x > of f(x). A first 
region, for a%x <C 1, characterized by the linear behavior 

f(x) ~ 1 - aix , (3.3) 



like the exponential exp(— ai x) does for x — > 0. A second intermediate region, 
for x <C 1 <C aix, where /(x) is characterized by the power law behavior 

/(x) ~ x" si , (3.4) 
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with slope s\ = l/(|^i| — J~i). 

Finally, for i> 1 we obtain the asymptotic power law behavior 



/(*) 



~ x 



(3.5) 



whose slope is now s 2 = l/d^al ~ r 2)- 

Thus, f(x) behaves like a power law function both in the middle and in the 
far region of x > with slopes Si and S2, respectively. In this sense, we call 
Eq. (13.21) a two-power law function. 

Form the above analysis, we easily realize that the constant ai, introduced in 
the definition of U ai (x), gives approximatively the width of the intermediate 
region having slope Si. 

As an example, let us specialize Eq. (13.21) to the case r = 0. In this situation, 
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Fig. 1. Log-log plot of the two-power law function (13. 2\\ (solid line). The dot-dashed 
line is the cumulative integral of the log-normal function. The dashed lines denote 
the asymptotic extension of the function P{x) in the two-power law regions. 



the generalized exponential and logarithm assume, respectively, the expression 
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exp {K} (x) 



KX + \/T + 



K 2 X 2 
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(3.6) 



and 



r 



— ft 



2k 



(3.7) 



In figure 1, we plot the function 



P(x) = exp {Ki} (ai ln {Kl} (exp {K2> (-x))) 



(3.8) 



for the values K\ = 0.35, k 2 = 1.4 and a\ = 10 3 . In the same graphic, the 
dot-dashed line depicts the cumulative integral I(x) = f£° p(y) dy of the log- 
normal distribution 



The dashed lines represent the asymptotic prolongation of the power law be- 
havior of P(x) whose slopes are given, respectively, by si = = 2.85 and 
S2 = I/&2 = 0.71. We observe a good agreement between the functions I(x) 
and P(x) only in the low region of x. 



4 Application to econophysics 

In the following, we employ the function derived in the previous Section 3 to 
fit some distributions data obtained in the economy framework. 
We pose P(x) = /(— /?|x| M ), with (3 and /i fitting parameters, the cumula- 
tive distribution representing the probability of finding a value X equal to, or 
greater than x. 

In figure 2, we present the results of the fit (in log- log scale) for the data of 
the inverse cumulative distribution of the personal income of Japan (1975) 
obtained in |18j and USA (2000) obtained in [3], as well as, the data of the in- 
verse cumulative distribution of the gross domestic production of Brazil (1996) 
and Germany (1998) obtained in [3]. 

In every graphic, we report the dashed lines representing the asymptotic be- 
havior of P(x) in the two power law regions with slope given by S\ = fi/k\ 
and s 2 = ^/k 2 , respectively. 
The data fit are reported in table 1. 




(3.9) 
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Fig. 2. Log-log plot of personal income distribution for Japan (1975) [18J and USA 
(2000) [3] and gross domestic production distribution for Brazil (1996) and Ger- 
many (1998) [3j. The solid line represents the fit obtained with the two-power law 
function (j3.8j) . The straight dashed lines are plotted for convenience to indicate the 
asymptotic power-law prolongation. 



Table 1. 

Parameters for the cumulative distribution P(x). 



Country 




^2 


a x 


V- 






Japan (1975) 


1.14 


2.00 


390 


3.00 


3.5 


■ 10~ 4 


UK (1998) 


1.70 


0.75 


8 


2.12 


2.5 


■ io~ 4 


Brazil (1996) 


2.20 


1.53 


2 ■ 10 4 


1.99 


2.3 


■ io- 8 


USA (2000) 


2.00 


0.65 


231 


1.44 


6.0 


■ 10~ 6 



The crossover between the first and the second power law region, causing a 
reduction of the slope, with s 2 < 8\ (UK, Brazil and USA), is named kink 
effect [3]. Similarly, the crossover between the first and the second power law 
region causing an increase of the slope, with s 2 > s% (Japan), is named ankle 
effect. 
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5 Generalization 



Let us briefly discuss the generalization of the method introduced in Section 
3 in order to generate functions with more than two power law behavior. This 
can be accomplished starting from the building block function 

n K} (x) = exp {K . r . } (ai \n {Kun} {x)) , (5.1) 
and introducing the quantity 




Personal Income 



Fig. 3. Log-log plot of the cumulative distribution of Japan for the year 1988. It is 
observed a deviation from the Pareto behavior in the highes income region (x > 10 4 ). 

where a = (01, 02, • • • , Cn-i) is a (n— l)-vector whose ith entry = (k^, r i: a^) 
contains the relevant informations about the slope and the width of the ith 
power law region. It is easy to verify that the function 

f{x) = LT M o exp {Kn rn} (-x) , (5.3) 
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exhibits a n-power law behavior. 

In figure 3, we report the fit of the 1998 Japanese income data obtained in 
[T9"j by employing the function 

P(i) = %oexp w (-/?/) , (5.4) 

derived from Eq. ( 15. 3ft for n = 3 and = 0. The fitting data are K\ = 
0.71, k 2 = 1.12, k 3 = 2.77, ai = 10, a 2 = 4 ■ 10 5 , = 6.00 ■ 10~ 9 and y. = 2.00. 



6 Conclusions 

We have derived a simple method which permits to generate functions with 
a multi-power law behavior starting the deformed logarithm hxj KjJ .\(x) and 
the deformed exponential exp| K r |(x), recently derived in [FT], which exhibit 
a single power law profile. An explicit two-power law function has been con- 
structed starting from the ^-exponential and its inverse, the ^-logarithm. We 
have employed this function to fit the inverse cumulative distribution of the 
personal income and of the gross domestic production of several countries, 
showing a good agreement among the analytical and the empirical data for a 
wide range of values. 
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